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Systems of one-dimensional magnetic or electronic chains coupled to so-called ladders have 
attracted increasing attention as a scenario in which two-dimensional physics can be studied by 
starting from well established one- dimensional theories fTj-^. The properties of the undoped 
systems, i.e. spin ladders, have been thoroughly investigated by a variety of methods |l|-§|] and 
recently integrable models with tunable interaction parameters have been introduced |7|-[T3|. To 



understand the properties of the doped ladder systems the dynamics of the holes created under 



doping in the background of antiferromagnetically interacting spins has to be studied [0 . This 
question has been studied within a model of coupled t-J chains forming a ladder-like structure 
in a series of papers |[T^-[16| with interesting results. Within this approach, however, multi hole 



scattering processes or the dynamics of anti-bonding hole states could not be accessed. 

Our aim here is to propose a modified t-J ladder model, which becomes integrable at certain 
limits and therefore allows us to investigate the properties of this system exactly and in more 
detail. In the thermodynamic limit, we can analyze further the nature of the ground state 
and predict the onset of excitations, thereby identifying the phase diagram of the system. The 
system we consider consists of two coupled t-J chains labeled a = 1, 2. The electrons with spin 
a on site j of the chains are described in terms of canonical Fermi operators c,° . In terms of 
these the chain Hamiltonian is HI'Jj = ^j{ Hl_j \ with 
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where V projects out double occupancies, S^" are spin operators on site j and rij = n^'t' +ti". 
is the total number of electrons on site j. 

The Hamiltonian of the ladder model we consider then reads 

H = J2 Hi- J + i^int + i^rung - f^N . (2) 

a 

Here, n is the chemical potential coupling to the total number of electrons in the system, i^int 
is the interaction between neighbouring rungs of the ladder, given by the expression H^^t = 
— "^A H\_j > < Hl_j \ , and if rung contains all the interactions present on a single rung 
of the ladder and reads 

ft™, = E [-tr {Y.fS'-^'S + "-) P + J' (s«sf - in«nf' ) + V'nfnf I . (3) 

This is similar to the terms in (|l|) describing the interaction of electrons on the bonds within a 
rung but contains independent coupling constants t', J' and an additional Coulomb interaction 



of strength V. Below we shall use a basis of states on a single rung which are eigenstates to @: 
at half filling (one electron per site) the electrons can form a spin singlet \s) = 4j(| ti)"! it)) 
or triplet state |t+) = | tt) stc. (|cTr) = Co- Cr |00)). Doping of the system creates the states 
\a±) = -j^ilcrO) ± I Oct)), with a =t or |, describing a single electron with spin a in a bonding 
(anti-bonding) state or the Fock vacuum (two-hole state) \d) = |00) on a rung. Below we 
shall be particularly interested in the strong coupling regime J' ^ 1, V ^ fi + \t'\ near half 
filling. This implies that the triplet states \ta) in the two-particle sector will be energetically 
unfavourable. 

To study the phase diagram of the system at low temperatures we note that by excluding 
the triplet states and choosing the coupling constants on the legs of the ladder as J = 2t = 2 
the Hamiltonian (^ can be rewritten as 

5 

j 1=1 

where Ni, I = 1,2 (3, 4) is the number of bonding (anti-bonding) single particle rung states with 
spin t, J, and N^ is the number of empty rungs (on the remaining Nq = L — ^^ A*"; rungs the 
two electrons form a singlet). The (graded) permutation operator 11^^ interchanges the states 
on rungs j and k giving a minus sign if both rungs are fermionic, i.e. singly occupied. In the 
restricted Hilbert space we are considering here we can absorb J' into the definition of V which 
allows to express the potentials Ai in terms of the coupling constants in (|^) as 

Ai=A2 = fi+=t'-fi + V, 

A3 = Ai = fi^ = -t'-fi + V, A^ = V = -2/i + V (5) 

The relative strength of these potentials will determine the ground state and spectrum of low 
lying excitations in the ladder. 

Note that the Hamiltonian (^) is that of a gl{2\A) invariant super spin chain in the presence 
of external symmetry breaking fields and the grading of the states |s) and \d) is bosonic while 
the two doublets \a±) carry fermionic grading. This model is integrable by means of the Bethe 



Ansatz (BA) (see e.g. p7|-p0|). In the strong coupling regime near half filling, where the 
ground state shows little deviations from the dimerized state ®j\s)j doping will lead to the 
condensation of either the fermionic single hole states |cr-|-) or the bosonic double hole state 
\d) into the ground state - depending on the relative strength of the chemical potentials /i_, 
/i+ and the rung interaction V . For simplicity we shall consider the case fi+ = fi- = fl below, 
i.e. two degenerate single particle bands. Due to the grading different sets of BA equations 



have to be considered to describe these processes (see e.g. [|T^,0,^): since the Bethe Ansatz 
states are highest weight in the gl{2\4:) algebra we have to deal separately with two regions 
in the low temperature phase diagram of the t-J ladder (^ between half- and quarter-filling, 
corresponding to 4A^5 < '^i^i Ni and AN^ > J2i=i ^h respectively. 

For the first case we have to order the rung basis as |s), \<J±), \d) and the eigenstates of (§) 
are parameterized by solutions of the BA equations corresponding to this grading (BFFFFB) 

r Ma 
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u^2{\^'-^') = n^i(4'^^-^r'')n^i(^r-^^i^^'^)' r=2,...,4, (6) 
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with e„(x) = (x + in/2)/{x — in/ 2). Here the number of roots on the rth level is related to the 
numbers Ni in (^) by M^+i = Mr — N^ (we define Mq = L). The energy of the corresponding 
state is 

Afi -. 5 

In the thermodynamic limit L —>■ oo with Mr/ L fixed the solutions of (|^) are built on two types 
of so-called strings, namely complexes 

^^^_„) ^ ^^r) ^ 1(^ _ 2^) ^ A; = 0, . . . ,n , < n < r and 2 < r < 4 

corresponding to bound states of holes with different spin and parity with densities p^*"^ (^) and 

Ai'') = A(;) + l(m + 1 - 2A:) , fc = 1, . . . , m , 2 < r < 4 

corresponding to bound spin states with densities am (A), r = 2,3,4. On level 1 and 5 all 
solutions are real with densities aj'' (A^*"^), r = 1,5. 

By minimizing the free energy of the system (see e.g. JIB|,pD|) one finds that at most the 
bound states of single particle rung states \a±) with density p*-^-' and empty rungs \d) with 
density a\ are present in the ground state. The densities are determined through the Bethe 
Ansatz integral equations 

\ ) ( ^\ —- n A ^\ — I Alt T^ ( ^ — ii\ rv- H ii\ — I Alt n -, ( ^ — it\rr^ ^ 



p^>{x) = ai{x) - dy K{x-y)p^ >{y)- dy ai{x-y)a\\y) , 

(8) 
.P(.) - faya.ix-ymy) 

J s 

3 



where K{x) = X]n=i '^2n(a^) and an{x) = 2n/ {ixi^Ax^ + ra^)). The corresponding dispersions e^ 
and ed are given by similar integral equations 

es{x) = 27Ta4{x) -A{2 + fi) - dy K{x -y)es{y) + dy ai{x -y)ed{y) , 

(9) 
ed{x) = ft - V - dy ai{x - y)e,{y) . 



The chemical potential and the rung interaction determine the ranges of integration 

dy =\ + \dy 

through the condition ea(±(5a) = 0. In terms of the solution of (|]) the density of holes is 

2Jdyp^'\y)-jdyaf\y). (10) 
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For V < jl < —2 both e^ and e^ are positive and hence the dimer state is the ground state 
of the system. For larger chemical potential single hole states \(j±) with dispersion e^ begin 
to condensate into this ground state. The electronic properties in this phase are those of the 
degenerate supersymmetric t-J model [^. For V > p, — f dy ai{y)es{y) the system enters a 



phase where in addition double hole rung states \d) with energy e^ will be present in the ground 
state. 

For ft <V we have to order the the rung states as \s), \d), \a±) basis leading to a different 
set of Bethe Ansatz equations corresponding to BBFFFF grading (Mg = 0). 

^ Ml M2 

k=/=j k=l 
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.(2) \(i)\ _ TT . m(2) 



1=1 m=l 

Mr Mr-1 Mr+1 

n^^(^r-^i^^) = n^i(4^'-^r'')n^i(^r-^^^'''^)- -=3,..., 5. 

kf^j 1=1 m=l 

(Note that due to the reordering of the basis the number of double hole rung states is now A^i.) 
The energy of the corresponding Bethe Ansatz eigenstates is again given by (|^) with A\ = V 
and A; = /i for 2 < / < 5. Similarly as for (^ the solutions of these equations can be classified 
in strings with densities Om for positive integers m and r = 1,3,4,5 and single hole bound 
states with densities p*-^-*, r = 3,4,5. Only real roots with density a\ are possible on the 



second level of (|TT]). The ground state configuration is completely determined by 



p(5)(x) = - dyK{x-y)p'^^\y)+ dya4{x-y)a'^\y), 



(12) 



a[^\x) = ai{x) + dy a^i^x - y)p^^\y) - I dy a2{x - y)cT^i^\y) ■ 

J s J d 

The dispersion of these excitations are determined through the integral equations 
e^(x) = 4:{V - p) - dy K{x-y)es{y)+ dy ai{x-y)ed{y) , 

J s J d 

ed{x) = 27rai{x) - {2 + V) + dy 04(0; - y)es{y) - dy a2{x - y)ed{y) . 



(13) 
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and the density of holes is rih = J^dy crj (y) — 2 J^ dy p'^^\y). 

With these equations the phase diagram of the t-J ladder @) can be completed: for ji < 

V < —2 the system is in the dimer phase, increasing the rung interaction V the system enters a 
phase where the ground state is a Fermi sea of double hole states with dispersion e^ moving in 
the background of dimer states \s). Increasing fl in this phase beyond | i7V — 2j both double 
and single hole rung states are present in the ground state configuration. 

Numerical integration of the integral equations for the densities and dressed energies allows 
to determine the phase boundaries as a function of the hole density n^ and the rung interaction 

V = 2jl — V in Eq. (^). In Fig. |I] the resulting phase diagram is shown. For V ^2 condensa- 
tion of the single hole states \a±) is energetically favourable. For sufficiently strong attractive 
rung interaction, however, doping of the system first introduces double hole rung states \d). 
Interestingly, we find that ioi 2 > V > 2 — ^{tjj (|) — ip (|)) ^ 1.676 (^(x) is the digamma 
function) increasing the hole concentration can drive the system from the latter phase into the 
one without double hole rung states. 

In summary we have studied the phase diagram of a system of coupled t-J models in the 
framework of an integrable system. By choosing properly the interaction between the neigh- 
bouring rungs the resulting ladder model can be brought into the form of a (7/(2|4)-invariant 
exchange model in the presence of symmetry breaking fields, which enter the Hamiltonian as ki- 
netic and interaction terms on a single rung. The spectrum of excitations and zero temperature 
phase diagram have been determined by means of the Bethe Ansatz. Extensions are possible 
by assigning different energies to the single particle rung states of different parity {t' 7^ 0) and 
including the triplet rung states (J' finite). As has been shown for the undoped case new phases 
can be expected for finite J' |T0|. 



HF acknowledges partial support by the Deutsche Forschungsgemeinschaft under Grant 
No. Fr 737/2 and AK acknowledges with thanks the Fellowship Grant from the Alexander von 
Humboldt Foundation. 

References 



[1] 

[2] 
[3] 
[4] 

[5] 

[6] 



[7; 

[8 
[9 

[lo: 

[11 
[12: 

[13 

[14 

[is: 

[16; 

[17 
[18 
[19 
[20 

[21 
[22 



E. Dagotto, J. Riera, and D. Scalapino, Phys. Rev. B 45, 5744 (1992). 

T. Barnes, E. Dagotto, J. Riera, and E. S. Swanson, Phys. Rev. B 47, 3196 (1993). 

S. Gopalan, T. M. Rice, and M. Sigrist, Phys. Rev. B 49, 8901 (1994), [cond-mat/9312026i 

S. Brehmer, H.-J. Mikeska, and U. Neugebauer, J. Phys. Condens. Matter 8, 7161 (1996), 
|cond-mat/9604136| . 



S. Brehmer, A. K. Kolezhuk, H.-J. Mikeska, and U. Neugebauer, J. Phys. Condens. Matter 
10, 1103 (1998), |53Hd'-mat/9710111i 



A. K. Kolezh uk and H.-J. Mikeska, Int. J. Mod. Phys. B 12, 2325 (1998), |5mld 
mat/9803176| . 



V. Yu. Popkov and A. A. Zvyagin, Phys. Lett. A175, 295 (1993). 

H. Frahm and C. Rodenbeck, J. Phys. A 30, 4467 (1997), |cond- mat/ 9 70^1153 



S. Albeverio, S.-M. Fei, and Y. Wang, Europhys. Lett. 47, 364 (1999), [coiid-mat/ 9807341 . 
Y. Wang, Phys. Rev. B 60, 9236 (1999), |cond-mat/9901168i 



N. Muramoto and M. Takahashi, J. Phys. Soc. Japan 68, 2098 (1999), [coiid-mat/ 9902007 
M. T. Batchelor and M. Maslen, preprint , |cond-mat/9907480| . 



M. T. Batchelor and M. Maslen, J. Phys. A 32, L377 (1999), |cond-mat/ 9907134 , 



S. Gayen and I. Bose, J. Phys. Condens. Matter 7, 5871 (1995). 

I. Bose and S. Gayen, J. Phys. Condens. Matter 6, L405 (1994). 

L Bose and S. Gayen, J. Phys. Condens. Matter 11, 6427 (1999), |cond-mat/981214q . 

P. P. Kulish, J. Sov. Math. 35, 2648 (1986), [Zap. Nauch. Semin. LOMI 145, 140 (1985)]. 

P. Schlottmann, J. Phys. Condens. Matter 4, 7565 (1992). 

F. H. L. EBler, V. E. Korepin, and K. Schoutens, Int. J. Mod. Phys. B 8, 3205 (1994). 

P. Schlottmann, Int. J. Mod. Phys. B 11, 355 (1997). 

F. H. L. EBler and V. E. Korepin, Phys. Rev. B 46, 9147 (1992). 

A. Foerster and M. Karowski, Nucl. Phys. B 396, 611 (1993). 



V 



4.0 



2.0 



0.0 



-2.0 



-4.0 





CJ> 




^^^ 






dimer 


d> 


"'x,,^ 






- 



0.0 



0.5 



n, 



Figure 1: Phase diagram of the t-J ladder with rung interaction @ for t' = and J' ^ 1 as 
a function of the hole concentration n/^: For large repulsive V the ground state is a Fermi sea 
of single hole states \a±) propagating in the background of rung dimers. For sufficiently strong 
attractive rung interaction the double hole rung states \d) are energetically favourable. In the 
intermediate region (shaded) dimers and both types of hole rung states coexist in the ground 
state configuration. 



